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Abstract 

The calculation of both resistive and ideal plasma equilibria amounts to solving a 
number of magnetic differential equations which are of the type B ■ V$ = s. We apply 
the necessary and sufficient criterion for the existence of the potential $ and find that 
a static equilibrium configuration of a magnetically confined plasma does not exist in 
axi-symmetric toroidal geometry. 
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I Introduction 

The term 'magnetic differential equation' was coined by Kruskal and Kulsrud [1] in 
1958. It arose from an attempt to model a 'resistive plasma equilibrium' with the 
equations: 

jxB = Vp (1) 
rotB = w (2) 
divB = (3) 
E + vxB = r]j (4) 
rotE = (5) 
From (|SJ) follows E = — V$ and from (£Q) after taking the scalar product with B: 

B -V$ = -r]j ■ B = s (6) 

This equation was called a magnetic differential equation. Newcomb [2] formulated a 
criterion for its solvability in 1959: 

Sdl (7) 



B 
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where the integration has to be carried out around closed field lines. 

Equations (1-5) are the basis for calculating the classical diffusive particle losses in 
a magnetically confined plasma. Pfirsch and Schliiter [3] used a magnetic model field 
which did not exactly satisfy equations (1 - 2), but they obtained, nevertheless, an 
estimate for the diffusion coefficient in a toroidal configuration. Later on Maschke [4] 
generalized their results for a true equilibrium configuration and confirmed the order of 
magnitude of the predicted losses. The present state of the art to calculate collisional 
transport in tokamaks including the effect of trapped particles can be found, e.g., in 
[5]. 

In this paper we analyze again the solvability of © and find that Newcomb's crite- 
rion (J2J) is necessary, but not sufficient. In general, the integral § V(f> ■ dl taken around 
any closed loop, which does not surround the central hole of the torus, must vanish, if 
the potential <fi is to exist. If we restrict ourselves to axi-symmetric configurations (e.g. 
tokamak), we see that (J7J) is only applied to a special class of closed loops, namely to 
the nested contours of the magnetic surfaces. In Section IV we apply the necessary 
and sufficient condition © on © and find that the inhomogeneous part s of the mag- 
netic differential equation © must vanish. Since s is given by other requirements, the 
condition for the existence of the potential cannot be satisfied. 

In Section V we reduce the force balance JI} to a magnetic differential equation 
which is of the same structure as ©• It turns out that the criterion for its solvability 
can also not be satisfied in general. The reason is that the transformation properties 
of the cross-product in are at variance with the transformation properties of the 
pressure gradient. As a consequence we come to the conclusion that the assumed 
stationary state of a magnetically confined plasma as described by (1 - 5) does not 
exist. 



II Equilibrium configuration 

We introduce a cylindrical coordinate system with unit vectors (e*R, e^, ez) and take 
the toroidal angle as an ignorable coordinate (axi-symmetry). The magnetic field 
satisfying (jSJ) may be written as: 

- F 

B = B p + By = Vif> x Vip + — e v (8) 

K 

where if) is the poloidal magnetic flux function. From (JI} and the poloidal component 
of © follows then: 

p = pW , F = F(ij) (9) 
and the current density may be expressed by 

F' - ( FF'\ 
j = —B p +[Rp' + — )e ip (10) 

where ' denotes differentiation with respect to if) . Inserting this into the toroidal 
component of (J2J) yields the Lust-Schluter-Grad-Rubin-Shafranov equation [6]: 

d 2 ip 1 dip d 2 ip / 2 / pe ,a nn 

dv-RdR + dz-i = AllJ = -y Rw+FF ) 

This equation is usually solved numerically with suitable boundary conditions to find 
the equilibrium magnetic field in axi-symmetric plasma configurations. There exist 



also analytical solutions (see Appendix A) which contain free parameters so that the 
solution can be adapted to realistic boundary conditions. Recently equation (|TT|) has 
been modified [7] to include approximately the effect of magnetic islands which break 
axi-symmetry. At this point, however, we assume that a solution of (fTTj) is given for a 
particular choice of p' and FF' . 



Ill Standard method for calculating the electric field and the velocity field 

The method used in [3] and [4] for calculating the electric field is straightforward. 
Because of (J5J) the electric field may be written as: 



E = E P + E V = -V0 + -— e v 

LTlti 
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where U is the loop voltage produced by a transformer. Substituting this into (@J and 
taking the scalar product with B yields the magnetic differential equation: 



B P ■ = ( V j v -E 9 )^+ ^B 2 p 



(13) 



with the solution: 



<p = O (i/j) 



lo 



f n F' 

- — Bl 



fJ-o 



Rdl 
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where dl denotes a line element on a contour of a magnetic surface if) = const. In order 
to obtain a single- valued potential, the closed integral around a magnetic surface must 
vanish: 



F ^ rjF' 



E„) - + ^—Bl 



Ik) 



Rdl 
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This equation reflects Newcomb's condition (JJJ) and puts a constraint on the choice of 
the functions p' and F' to be used in (fTTj) . 

The velocity field is obtained by taking the cross-product of (jlj) with B: 



v = a 



(R, Z)B+(E-nfj x B/B 2 



(16) 



The arbitrary function a may be expressed by the divergence of the velocity, thus 
resulting in a second magnetic differential equation because of (J3*|): 



B p -Va = div (v- (e - rjj) x B / B' 



(17) 



The solution is: 



a = (Xq 



(V>) + / divfv- (E-nj) x B / B 2 ) 



Rdl 



with the integrability condition: 



<p div ■ 



,Rdl 



div ([E -rjj) x B/B 



>2> Rdl 
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which is usually expressed in the form: 

Rdl , , f ( => - „ \ „ , Rdl 



div v 



(20) 



by application of Gauss' theorem. 

It appears now that the problem is solved, first by prescribing a flux function 
satisfying (fTTj) and (|15p. a divergence of the velocity subject to the constraint (f2T)l) . 
and secondly, by evaluating (fH| an d (fTS|) . The solutions may be substituted into 
(fT2*|) and (fTBj) to yield the vector fields E and i7. The integration functions <fto and 
Q!o describe an arbitrary rotation in poloidal and toroidal direction. In the following 
it is shown that this expectation is not justified, since (|15|) and (J2U|) turn out to be 
necessary, but not sufficient conditions. 



IV Application of the necessary and sufficient criterion for the existence of 
the potential 

Before we apply the necessary and sufficient condition (jSJ) for the existence of the 
potential <f), we write f!13|) in the form: 



Er B r + E z B z 



F n F' 
(VJv - E v ) - + V —Bl 



(21) 



and take the gradient of this equation: 



„ dE R n dE z „ dB R „ dB z ds 

td ££ i r> ±_ i jp tz. i Tp a 
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dR dR 
dB z _ ds_ 

" dZ 



dZ * dZ n 8Z " dZ 
Together with (j2~TJl . (jUJ), and the necessary and sufficient condition (J5J): 

dE R dE z 



(22) 



(23) 



one obtains from 
nents: 



dZ dR 

two magnetic differential equations for the electric field compo- 



B p -V[^- 
p \RB 2 



LJL(JL\ b v( Ez 

RdR \B Z ) ' P \RB R 



(24) 



In analogy to ([H)) they have the solutions: 



RB* 



Ez = R Br 
Inserting this into (}2*T]) yields: 



d 



dz \b r J |w| 
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Since the lower bound Iq of the integral may be chosen arbitrarily close to the upper 
bound I, the integral can be made to vanish so that equation (|25Jl leads to the condition: 

hM = -fiM = <t>oW , 5 = (27) 

as the right-hand-side of (J2H)) is not a function of ijj alone. Result (|2~7|) can also be de- 
duced from the transformation properties of the electric field vector whose components 
are given in (|25[). This is demonstrated in Appendix B. 

The local condition (|27|) is more restrictive than Newcombs integral condition (|15|). 
and it is not compatible with s as given in (|2*T|). Hence, we must conclude that the 
potential does not exist. 



V Magnetic differential equations and the force balance 

In case of the static force balance (0) scalar multiplication with the magnetic field 
yields the homogeneous magnetic differential equation: 

B ■ Vp = (28) 

which satisfies condition ()27|) so that an ideal equilibrium with rj = should be possi- 
ble. The necessary and sufficient condition for the existence of the pressure 'potential' 
is, however, that the curl of (JTJ) vanishes: [p • V) j = (j ■ V) B, which again leads 
to an inhomogeneous magnetic differential equation for the toroidal component of the 
current density together with (J2J): 

As shown in the previous Section, it has only the solution j v = Rf (i/j) so that the 
radial component of the current density must vanish. 

This may also be shown by considering the so called 'force-free' situation where the 
pressure gradient vanishes and ([lip becomes: 

A*i{; = -FF' = -g (ip) (30) 

Applying Stokes's theorem on this equation by integrating the toroidal current density 
over the area enclosed by a magnetic surface one has: 

^ [W^ „ If , ,s dRdZ 

With: 



one obtains from (j30p : 



, R di = g g ^)^- (si) 

1>= [ dF , W = - VF (32) 

J 9(F) 9 



A*F = -^-^\VF\ 2 ^-[-\ (33) 



F F 1 1 dF 
Stoke's theorem applied on this equation gives: 



Substitution of into (j3*Tj) yields on the other hand: 



F r \VF\ „ £/ di2 

dl=Qg —et- (35) 



g J R jj R 

Elimination of the line integral over the poloidal current density on the left-hand sides 
of an d (|3*3j) . and using (|3*2|) again results in an integral equation: 

rr dRdZ ^ U (g 2 1 / o ^ dg\ t ^ Il2 \ dRdZ . . 

which can only be satisfied for g = FF' = in agreement with (J27|) and (|29jl. We 
demonstrate this explicitly in Appendix A by adopting a 'Soloviev solution' [8] with 
g = const. 

In view of this result it becomes doubtful whether the condition p = p {ip) resulting 
from (J28j) in axi-symmetry defines a scalar pressure, or, in other words, whether the 
cross-product j x B, which is in principle an antisymmetric second-rank tensor, can 
have the same transformation properties as the polar vector field Vp, at least under 
certain circumstances. 

In order to investigate this question we apply the Laplace operator in Cartesian 
coordinates on 



B ■ V (Ap) + Vp • AB 

^ 2 / dB x d 2 p ^ dB y d 2 p ^ dB z d 2 p ^ ^ d 2 p ^ ^ d 2 p ^ ^ d 2 p \ _ q /gy\ 
\ dx dx 2 dy dy 2 dz dz 2 dxdy dxdz dydz J 

. dB x dB y dB x dB z dB y dB z 

kl = ~7T~ + a ' k2 = ~1T~ + ~7T~ ' k * = -7T- + -Tr- 
oy ox dz ox dz dy 

The gradient of 

Ox \ ox Oy Oz J 

8/ 4 + 4 + 4)^0 (38, 



dy \ dx dy dz 

A ( Bx ^P +B ^P + Bz ^l 
dz \ x dx v dy z dz 



yields three equations which can be used to eliminate the mixed derivatives in ([37)1 . 
The bracket may then be written in the form: 



/ fc 3 .R, h\ g ( k 2 B y _ h \ ^ f k x B z _ h_ _ 

\B y B z By B z J \B X B Z B x B z ) y \B x B y B x B y j 



(39) 



It turns out that this expression does not have the transformation properties of a 
scalar, in contrast to the first two terms in (|37j) . By 'scalar' we refer to a quantity 



which does not change its value, when it is expressed in different coordinate systems 
as a function of space: p(x) = p(x'), x' — x, x\ = a^x^- The reason for the 
'non-scalar' property of is that the directed quantity I cannot be considered as 
a vector field which maintains its modulus as a scalar, when it is transformed into a 
rotated coordinate system. The inner product of the last term in (J39j) will, therefore, 
depend on the orientation of the coordinate system which would not be the case for 
the invariant inner product Vp • \B ■ B, e.g., or for the first two terms in ([37)1 . 
Similar remarks apply to the first three terms in (}3"9"j) which resemble the Laplacian 
of the pressure, but the second derivatives have coefficients which are all different so 
that the sum of these terms is not invariant, as compared to the Laplacian of a scalar 
field. Consequently, equation leads to an incongruity, when it is transformed into 
a rotated coordinate system, as it will depend explicitly on the rotational angle. 

In order to show this we choose a coordinate system which is rotated around the 
y-axis by an angle a: 

x = x' cos a — z' sin a , z = x' sin a + z cos a , y = y' (40) 

The transformation rules are: 

d d . d d . d d d d , _ N 

— = cos a — - sin a— , — = sin a— + cos a— , — = — (41) 
ox ox oz oz ox oz oy oy 

Applying these to ()39j) one finds that the transformed expression contains not only the 
components of the magnetic field and of the pressure gradient in the primed system, but 
in addition the rotational angle a. This may be most readily verified by transforming 
the coefficient l y . It should be invariant, since the ^/-component of the pressure gradient 
does not change under the assumed rotation so that the second term in (}3"9"j) should 
be the same in the rotated system. Instead one obtains: 

, _ (k' 2 B y , - k[B z , - k' 3 B x ,) cos 2a + (k[B x , - k' 3 B z , + (gf - B y ) sin 2a 
v (B x > cos a — B z i sin a) (B x i sin a + B z > cos a) 

(42) 

where the k' n are defined as in (|H7jl in terms of primed derivatives of the primed field 
components. For a = one returns to the expression l y as given in ()39|). Evidently, 
the expression (|39j) does not transform like a scalar field which is only the case when 
the pressure itself, as defined by (J2EJ) in the form p (<0) , is not a scalar, contrary to 
our assumption. In Appendix A we show this explicitly for the pressure as given by a 
Soloviev solution. 

From the transformation properties of (|39jl . which are a consequence of the hypo- 
thetical equation (JTJ), we infer that a static equilibrium configuration does not exist, 
since the transformation properties of the vector cross-product are incompatible with 
those of the gradient of a scalar pressure. 



Conclusion 

The stationary equilibrium equations (1 - 5) of a magnetically confined plasma may 
be formulated in terms of magnetic differential equations. We have shown that these 
equations have ambiguous solutions in axi-symmetry, unless their inhomogeneous part 
vanishes. Furthermore, due to the transformation properties of the cross-product j x B, 



it cannot be set equal to the gradient of a scalar pressure. As a consequence, the set 
of equations (1 - 5) is not solvable. 

The time dependent terms which are omitted in the stationary magneto-hydro- 
dynamic model can apparently not be neglected. Consequently, magnetic confinement 
of a plasma in toroidal geometry leads inevitably to temporal changes of the pressure 
and the electromagnetic field. This may not be necessarily a unidirectional temporal 
evolution, but turbulent fluctuations could lead to an average 'quasi-stationary' state, 
which would not be strictly axi-symmetric any longer. 
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Appendix A 

Soloviev has constructed a flux function [81: 



* = \ {coR 2 + bRt)z> + ^ ( R 2 - Rlf (A.l) 



2 

which substituted into (JSJ) yields the poloidal magnetic field components: 

% ( _ n2 , 7. r>2\ d „ v2 i a ~ C / D 2 D 2 



B R = -- (c R 2 + bR 2 ) , B z = c Q Z 2 + — ^ [R 2 - R 2 ) (A.2) 

When these are inserted into the toroidal component of (J2J) one obtains: 

b R 2 

(t j v = -aR--^- (A.3) 

so that (jlUj) and ([lip are satisfied with: 

p' = -H a, FF' = -bR 2 (A A) 

Recently, expressions similar to (|A.1|) were published in [9], which have more adjustable 
parameters. They allow to model tokamak plasma configurations more realistically 
than (|A~TJ) . 

With g = —bRl equation (JHtjj) becomes: 

and from (|A.4J) follows: 

F = ^F 2 ~2bR 2 ^ (A.6) 

Converting the left-hand-side of (|A.5|) into a line integral and performing a partial 
integration on the first term of the right-hand-side of (|A.5|) yields with (|A.6|) : 

R 2 



Ri 



R2 Z(R^) dR+ u(JJ%Z_d± + W\ 2 \ dRdZ 



Ri 



FR N \ F 3 dZ F 3 R 



(A.7) 



where R\ and R2 are the points where a magnetic surface cuts the mid-plane Z = 0. 
The Z - coordinate on a magnetic surface is expressed as a function of R and ip with 
(lA.ljl . Collecting terms equation (jA.7|) becomes with (JHJ): 

= Fg>ff(bRtZB R + R(Bl + Bl)) ^ (A.8) 

Inserting the magnetic field components as given in ()A.2|) with a = one finds that the 
double integral over the cross-section of the plasma inside a magnetic surface does not 
vanish which requires then g — to satisfy (jA.8|) . This result was already expected 
from the condition (J2TJ) and the magnetic differential equation (J29J) . 

Soloviev's solution may help to understand the conclusion following from (J42J) that 
B • = does not define a scalar pressure. When the first expression in flA.4|) is 
integrated, one obtains for the pressure: 

p = p -^ aijj (A. 9) 

Because of B = rot A the flux function ijj in (JHJ) is related to the toroidal component 
of the vector potential: ip = RA V . This expression may be considered as the Z - 
component of the vector field Ren x A^e^. An arbitrary function of a single vector 
component does, however, not transform like a scalar field. Writing equation (jA.9|) in 
the form: 

p = p - ii G aRA v = p - n Q a(xA y - yA x ) (A.10) 

and transforming this expression into a coordinate system which is rotated around the 
y - axis by an angle a as in Section V: 

x = x' cos a — z' sin a , z = x' sin a + z' cos a , y = y' 

(A.ll) 

A x = A x i cos a — A z i sin a , A z = A x > sin a + A z i cos a , A y — A y > 
one obtains: 

p = po — hq a [(x'Ayi — y'A z i) cos a + (y'A z i — z A y >) sin a] (A. 12) 

This expression contains not only the coordinates and the components of the vector 
potential in the primed system, but in addition the rotational angle a. Hence, the 
pressure obtained from (jA.9|) does not transform like a scalar. 



Appendix B 

The derivation of the condition s = in Section IV rested on a particular choice of 
the lower bound of the integral (J26)) . More generally, result (J27j) follows also from the 
vector character of the electric field whose components are given in (|25|) in cylindrical 
coordinates. If one formulates equation (|21j) in spherical coordinates (R = r sin 9 , Z = 

r cos 9 , f = ip): 

E r B r + E e B e = s (B.l) 

the method applied in Section IV yields now expressions for the field components in 
spherical coordinates. Taking the gradient of (|B.1|) one obtains with (JHJ) and (JHJ) two 
magnetic differential equations for the electric field components: 



which have the formal solutions: 



E r = R Bq 



Eg = R B r . 



1 d (sr\ dl 



1 d ( s 



rd6 \B r J |W| 



dl 



+ h W 



(B.3) 



Writing Eg = (s — E r B r ) j Bg and Ez — (s — Er Br) j Bz one obtains for the modulus 
of the electric field from (IB. 3D: 



E; • E 2 R 2 B 2 (f 3 I hf I 



and from we have: 



s 2 2sRB r . 



Bg 



(/a + h),h = J 



1 d ( s r \ dl 



lo 



r dr \B e J \Vip\ 



(B.4) 



E R +E 2 =R 2 B 2 p (f l + h)^ 



n s 2 2sRB R 



B 2 Z 



B, 



r// 



(B.5) 



Subtracting both equations one finds: 
(fi + hf - (/ s + I 3 ) 



f? 2 5 2 I £ 2 fi 2 



(B.6) 



If one differentiates this equation along the poloidal magnetic field lines, one obtains 
a linear relationship between the integrals 1\ and I3: 



h {h + h) = h (/1 + A) + k 2 



(B.7) 



h = B p - V/i - -B„ • V 



h = B p - V/ 3 - S p • V 



^5 e 5 2 



, 5 P • V/i 



]__d_ 

RdR \B Z 
1 d f ST 



Rr dr \B. 



k 2 = B p -V 



1 1 



2R 2 B 2 \B 2 Bl 



S f Br - B r - 

RWM B '- Vh -B e B >- Vh 



Eliminating from ()B.6|) and (jB.7|) the integral ^3 one obtains a quadratic equation for 
(fi + h): 



(./': 



A; 2 



2 + 



s / Br B r k± 



2sB r k 2 k\ 



RBj\B z B e k 3 
s 2 ( 1 \\ 



RBgBlk* k\ R 2 Bl \B 2 Bl 







(B.8) 



■p '"o '-j ~p 

Because of f|2K|) : fi + h = Er/(R Bz) equation (|B.8|) yields an explicit algebraic 
expression for Er as a function of s and the poloidal magnetic field components. This 



result is not compatible with Er given as an integral in the first equation of (|25|). 
unless condition (|27jl is satisfied. In this case follows from ()B.6|) : fi=fi = <fib 2 and 
the potential becomes a function of ifi only. 



References 

[1] M. D. Kruskal, R. M. Kulsrud, The Physics of Fluids 1 (1958) 265 
[2] W. A. Newcomb, The Physics of Fluids 2 (1959) 362 

[3] D. Pfirsch, A. Schluter, Max-Planck-Institut fur Physik und Astrophysik, Report 
MPI/Pa/7/62 (1962) (unpublished) 

[4] E. K. Maschke, Plasma Physics 13 (1971) 905 

[5] D. A. Gates, H. E. Mynick, R. B. White, Physics of Plasmas 11 (2004) L45 

[6] R. Lust, A. Schluter, Zeitschrift fur Naturforschung 12 (1957) 850 

V. D. Shafranov, Sov. Phys. JETP 6 (1958) 545; Zh.Eksp.Teor. Fiz. 33 (1957) 
710 

H. Grad, H. Rubin, Proc. 2nd U. N. Int. Conf. on the Peaceful Uses of Atomic 
Energy Geneva 1958, Vol. 31, 190, Columbia University Press, New York (1959) 

[7] X. Liu, J. D. Callen, C. G. Hegna, Physics of Plasmas 11 (2004) 4824L 

[8] L. S. Soloviev, Sov. Phys. JETP 26 (1968) 400; Zh. Eksp. Teor. Fiz. 53 (1967) 
626 

[9] C. V. Atanasiu, S. Giinter, K. Lackner, I. G. Miron, Physics of Plasmas 11 (2004) 
3510 



